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1. Introduction 

An interacting particle system may be informally described as 'a Markov process 
consisting of countably many pure-jump processes that interact by modifying 
each other's transition rates' ([24], p. 641). It is often assumed, as in [24] and 
most of [40], that these pure-jump processes live in a finite or at least compact 
state space; however, in applications this is not always the case. In the present 
work (in Section 2) we give a general construction for interacting particle sys- 
tems where we assume only that the constituent pure-jump processes live in a 
separable complete metric space that is not required to be compact. We allow 
jumps occurring at a given site to affect the states (not just the jump rates) at 
neighbouring sites. 

As well as a general existence result, we give general results on spatial limit 
theory. Typically, the interacting jump processes are indexed by the integer 
lattice with local interactions. This infinite system may be viewed as a limit 
of a sequence of finite systems, indexed by large finite 'windows' W„ C 7L d . 
Averaging a (time-dependent) quantity of interest over the finite system W n , 
one might expect to see it satisfy a law of large numbers and (functional) central 
limit theorem as W n t and we give general results of this type in Section 3. 

lattice and continuum models of sequential particle deposition. These have 
been much studied in the physical sciences literature [5; 21; 53; 61], and more 
recently in the mathematical literature [2; 6; 45; 46; 47; 49; 50; 51; 57]. In these, 
particles fall at random from above towards a surface (or substrate) and deposit 
themselves onto the surface (or fail to do so), in a manner which depends on the 
existing configuration of particles on the surface. According to the details of the 
model, the aggregation of deposited particles may be monolayer or multilayer, 
and the locations of deposited particles may be restricted to a lattice or may be 
continuous variables. 

The results of the present paper extend known existence and limit results 
to functional central limit theorems, and to deposition models allowing for dis- 
placement whereby an incoming particle may influence the positions of existing 
particles. Models of 'random loose sphere packing' [9; 10; 14; 55] often feature 
displacement in some form. The framework given here encompasses both lat- 
tice and continuum models, both in monolayer and multilayer form. Our limit 
results apply to quantities such as the total number of deposited particles, or 
(in the case of multi-layer models) the total height of 'exposed' particles (i.e., 
those which lie at the surface of the accumulation of particles) , and the average 
number of inter-particle contacts per unit area. Often in the physical sciences 
literature, such quantities are of interest, and their limiting values are assumed 
to exist without rigorous mathematical proof. 
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More generally, our results apply to any locally finite marked point process 
evolving jumpwise with locally determined jump rates and with locally bounded 
jump rates (here the word 'jump' could refer to births and deaths as well as 
movements of points) . We refer to these processes as spatial birth, death, migra- 
tion and displacement processes. They are described in general terms in Section 
4.1 and special cases such as deposition models are discussed in the rest of Sec- 
tion 4. In this type of continuum model, the 'state' at a site in 7L d refers to the 
configuration of points in a patch of M d . Our results add to previous work in 
[26; 27; 33; 52] on the theory of evolving point process. 

In recent complementary work, Garcia and Kurtz [25] define spatial birth- 
death processes in the continuum as solutions of so-called stochastic equations, 
and give criteria for existence and uniqueness; Qi [5 I] obtains functional central 
limit theorems for such processes (the work in [54] was carried out independently 
of that in this paper, which in turn was carried out independently of [25]). Unlike 
the present paper, [25] and [54] do not consder migration and displacement; on 
the other hand, their approach allows some relaxation of the assumption in the 
present paper of bounded range interactions. 

We discuss some discrete-space examples and open problems in Section 5. 
Subsequent sections are devoted to proofs. 

2. A general existence theorem 

In this section we give our general result on existence and uniqueness of inter- 
acting particle systems. Before describing this, we briefly review some concepts 
from the theory of time-homogeneous Markov processes. Suppose (E, £) is a 
measurable space. A probability kernel on E is a mapping \i : E x £ — > R such 
that n(-,A) is measurable for each A 6 £ and fJ>(x, •) is a probability mea- 
sure on E for each x £ E. A Markovian family of transition distributions on 
E is a collection of probability kernels (/i{,t > 0) with fJ,o(x, •) = S x for each 
x G E and with J E nt(y, A)/j, s (x,dy) = jji s +t{x, A) for all A € £, x, y £ E, and 
s,t > 0. The induced (transition) semigroup of operators (Pt,t > 0) is given by 
Ptf(x) = J E f{y)nt(x,dy), defined for all bounded measurable / : E — > K. The 
generator of this semigroup is given by 

Gf = limt- 1 (P t f-f), (2.1) 

defined for all bounded measurable / : E i— > M for which the limit exists with 
uniform convergence. Note that in this paper we require uniform rather than 
pointwise convergence here for Gf to be defined, as in e.g. [22] page 8 or [20] 
page 22. 

Suppose that E is a topological space and £ is the Borcl er-field on E. Suppose 
(fJ>t,t > 0) is a Markovian family of transition distributions on E, and (Qt,t > 0) 
is a filtration. A Markov family of processes in state-space E with transition 
distributions {fi t , t > 0) is a collection of stochastic processes > 0), defined 
for each x S E, which have right-continuous sample paths with left limits, are 
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adapted to the filtration, and satisfy P[£o = x] = 1 and P[^s+t G A\Q S ] = 
/j, t (^,A) almost surely, for each (x,A) € E x £. 

2. 1 . The general framework 

Let X be a Polish space (i.e., a separable, completely metrizable topological 
space). We consider Markov processes on the state space X v , where V is finite 
or countable (usually 7L d or finite subset thereof), endowed with the product 
topology. We assume that V is the vertex set of a graph V with 

D : = max{degree(u) : v £ V} < oo. (2-2) 

If the state of the process at time t is denoted £t, then £t(v) is the state at site 
v and lies in the space X; we shall refer to X as the local state space and to 
X v as the global state space. We do not require X to be finite or even compact; 
this feature puts us outside the scope of most of the existing interacting particle 
systems literature (though not all, as we discuss later). Many of the examples 
we consider are continuum systems where X is required to be infinite. 

Locally, the process evolves as a pure jump type Markov process. We assume 
finite range. That is, if a change occurs at site v £ V, then the instantaneous 
effect of this change is restricted to sites in a neighbourhood of v, and the jump 
rate at v is determined by the state of the system in a neighbourhood of v. 
The neighbourhoods are determined by our assumed graph structure on V; for 
each v £ V, define Af v C V to consist of v along with all adjacent vertices in 
the graph V. The neighborhood M v represents the set of all sites that can be 
instantaneously affected by changes at v or which affect the jump rate at site v. 

For each v £ V, assume we are given an operator G* which is the generator 
of a continuous-time pure jump type Markov process on the state space X^ 
with bounded rate function (for definitions see [38] page 238, or [24] page 635). 
Then for bounded measurable / : X^ v — > E, and for x € X^ , we have (see 
[38], Proposition 19.2, or [24], eqn 31.11 or [22], page 376) 

G*J(x)= I (f(y)-f(x))a* v (x,dy) (2.3) 

where the measure a*(x, •) is the jump rate kernel for the Markov process on 
X^ v (the generator G* generates the process as it would proceed if changes 
took place only 'at v'). Let a*(x) denote the total measure of a*(x, •), so that 
the next jump of the process with generator G*, when in state x, occurs at rate 
al(x) and is governed by the measure a*(x, -)/a^(x) when it does occur. We 
assume throughout that rate functions are uniformly bounded, i.e. that 

Cmax := sup{<(a;) :v€V,x€ X^"} < oo. (2.4) 

Let C be the class of bounded measurable functions / : X v — * M which 
depend only on finitely many coordinates. Let v £ V. We define a generator G v 
of a jump process on X v , as follows. Given x £ X v , and y £ X^ v , let x\v\y 
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be the element of X v which agrees with x outside Af v , and agrees with y inside 
M v . In other words, for w G V set 



(x\v\y)(w) 



y(w) if w S M v 
x(w) otherwise. 



For / G C, let the function f% : X M - -> R be given by 

fUy) := f(x\v\y). (2.5) 

Let x\n v be the restriction of x to J\f v (an element of X^ v ), and let ar|v\// v be 
the restriction of x to V \ Af v . Set 

a(/)(x) = G:(/^(xkJ). (2.6) 

Then G„ is the generator of a jump process with jump kernel denoted a v , given 

by 

a v (x,x\v\dy) = a* v (x\^ v , dy); 
a v (x,dz) = if z\ V \m v ^ x\v\AT v - 

Here, for z G X y the notation z\y\jf v ^ x\v\m~ v means there exists w G V \M V 
with z(w) 7^ x(iy). 

Given W CV, define the operator G 1 ^ on C by 

G w = J2 G *- ( 2J ) 

We are interested only in cases where W is finite or W = V. In the latter case, 
write simply G for G v . 

If W C V is finite, it is clear that G w is the generator of a jump process 
with jump kernel given by the sum over v G W of the jump kernels a v . Write 
for the associated semigroup. 



2.2. Existence and uniqueness 

Let the space X v be endowed with the product topology and Borel tr-field. Our 
first general result gives existence and uniqueness for a Markov transition semi- 
group on X v with generator G, and existence of a Markov family of processes 
in state space X v corresponding to this transition semigroup. This process has 
values in -D([0, oo), X v ) 1 and is called an interacting particle system, at least 
when X is finite (see [24], Chapter 32). Given a sequence of sets (W m ,m > 1) 
we write liminf m _ >00 (W / m ) for U~ =1 n£° =m W n . 

Theorem 2.1. (Existence and uniqueness result) Let G be the generator defined 
by (2.7) with W = V. Assume (2.2) and (2.4) hold. Then there exists a unique 
Markovian family of transition distributions (u t) t > 0) on X v , such that the 
associated semigroup (Pt,t > 0) has generator G on C. 
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Moreover, for any sequence (W m ) m >i of finite subsets ofV, and any t>0 
P t f{x) = lim P™ m f{x), V / eC,x e if liminf(W ro ) = V. (2.8) 

m — >oc m — >oo 

There also exists a filtration (Gt^t > 0) and a Markov family of processes 
G -X"^), adapted to (Qt,t > 0), wii/i transition semigroup (Pt,t > 0). 

Most previous existence results for interacting particle systems have been 
given for cases when the space X is taken to be finite (for example [24; 39]), or 
at least compact ([32], and Chapter 1 of [40]). Unlike in these results, we do not 
assume X is compact. If X is not compact then X v is not even locally com- 
pact, which also renders inapplicable much of the general literature on Markov 
process theory. Compactness of X would imply existence of at least one equi- 
librium distribution (see, e.g., [24]), but we do not here concern ourselves with 
equilibrium distributions. 

There is some previous literature in the case of non-compact X . In Chapter 9 
of Liggett [40], the space X is [0, oo) and the interactions have a linear structure. 
Our setup allows general noncompact X and does not require this linear struc- 
ture, although conversely, Chapter 9 of [40] does not require the interactions 
have the local structure considered here. Basis [7; 8] gives a general existence 
result in non-compact spaces, and Chen ([12], Chapter 13) gives a similar result, 
and also a uniqueness result. However, the general approach, and the detailed 
set of conditions, in those works seem very different from here. In particular, 
both [7] and [12] require a 'ground state' in X v to be specified, and take each 
approximating jump process, with semigroup (P t n ,t > 0), to have its state 
outside W„ to be given by the ground state (which can be viewed as a 'bound- 
ary condition'). With our approach, there is no need to specify a ground state 
(and the approximating jump processes are not needed). The examples consid- 
ered here are mainly in continuous spaces while those considered in [7; 12] are 
mainly in discrete spaces. 

The literature includes a number of papers on existence and other properties 
of particular interacting particle systems with non-compact X and non-linear 
structure (in some cases with unbounded jump rates). For example, [41] gives a 
construction for several such processes and the method can be adapted to others 
with non-linear rates (see [1]). In these examples the state space is not the whole 
of X v but a space of functions in X v satisfying a summability condition. Other 
examples include [58] and [4]. 

Our proof for Theorem 2.1 loosely follows the graphical representation of 
Harris [31] for finite-range particle systems, even though that paper, and subse- 
quent presentations such as those in [19; 28; 40], are also restricted to the case 
where X is finite. In extending this argument to a more general class of sets 
X, we use a measure-theoretic result from Kallcnberg [38] which is used there 
to express Markov chains as random dynamical systems. Using this result, any 
pure jump continuous-time Markov process with bounded jump rates in a Polish 
state space can be generated by a Poisson process in (0, oo) x [0, 1]. In proving 
Theorem 2.1, we shall show that the processes £f generated by G can be realized 
in terms of a family of independent Poisson processes, indexed by v 6 V. 
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3. Spatial limit theorems 

3.1. Further assumptions and definitions 

We use the following notation throughout. Given d £ N we let denote the origin 
in R d . For A cR d and u £ R d we write u + A for the set {u + u' : u' £ A}. For 
r > 0, let B r be the closed Euclidean ball of radius r centred at the origin in 
R d . 

For the results in this section, we make some further assumptions, in addition 
to those of Section 2. We assume that V = 1 d , that some specified 'window' Af 
around the origin (i.e. a finite symmetric set AT C 1i d with £ Af) is given, and 
that Af v is the set v +AT for each v £ Z d . We also assume translation invariance 
of the jump rates associated with each lattice point. In other words, we assume 
we are given a jump rate kernel a* (x, dy) on X^ , and that for v £ Z d , the jump 
kernel a* on X M - is the shifted version of a* given by 

a* v {x;T)=a*{L v (x),L v (T)), x £ X^" , T Borcl in X^, (3.1) 

where L v : X^ 1 ' — > X 1 ^ is given by 

(L v x)(w) — x(w + v), xeX^,weAf. 

We assume the total measure of a*(x,dy) is bounded, uniformly over x £ X^; 
this is enough to ensure that the condition (2.4) holds here. 

Let B denote the collection of all non-empty finite subsets of 1 . Let v be 
a probability measure on X. For A £ £>, let ,t > 0) be a realization of 
the X z -valued pure jump Markov process with generator G A given by (2.7) 
and with initial value = (£o'"(v),v £ 7L d ) whose entries ^' v (v) (v £ 7L d ) 
are independent and identically distributed states with common distribution v. 
Likewise, let , t > 0) be a realization of the X z -valued Markov process with 
generator G and satisfying (2.8) (as given in Theorem 2.1) with initial value £g 
whose entries £q (i;) (i; £ Z d ) are independent and identically distributed states 
with common distribution v. 

Our limit theorems are given in terms of stationary additive set functions 
on X z , which we define as follows. Given a measurable real-valued function 
H : X 1 ^ — > K, the corresponding stationary additive set function S^, defined 
for each A £ B, is a real-valued functional on the global state space X zd obtained 
by summing H over the states at (neighbourhoods of) sites in A, i.e., defined 
for x £ X v by 

sjHO = (*kJ)- (3-2) 

We shall give limit theorems, in the form of a law of large numbers (LLN) and 
a central limit theorem (CLT) as A becomes large, for the random variable 
Sff(£ t '"). The idea here is that the evolution of the interacting particle sys- 
tem is restricted to the finite set A, and the function H is also summed over 
(neighbourhoods of) states at sites in A. See Section 4 for some examples. 
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For A C Z d we write A c for Z d \ A, and define the ' neighbourhood', 'exterior 
boundary', and 'interior' of A, respectively, by 

Na ~ U veA Af v ; d oxt A := N A \ A; A° := 7L d \ N A «. 

Also let \A\ denote the number of elements of A. 

3.2. Limit theorems 

For 7 > and r > 0, we shall consider H and v satisfying the moments condition 
supE[\H(^ v W)n <oo. (3.3) 

We consider B- valued sequences (A n ,n > 1) satisfying the conditions 

\d eKt A n \/\A n \ -v asn^oo; (3.4) 
(Af An )° = A n Vn. (3.5) 

The conditions (3.4) and (3.5) hold, for example, if A„ = ([— n, n] n 2) d . 

We have the following LLN for functionals on interacting particle systems 
with finite range interactions. 

Theorem 3.1. (General LLN) Let d, M and H be as in Section 3.1. Let r > 0. 

Let v be a probability measure on X , and suppose that the moments condition 
(3.3) holds for some 7 > 1. Then for any B-valued sequence (A n ,n > 0) satis- 
fying (3.4) and (3.5), we have as n — > 00 that 

l^l" 1 ^ ^ EH(CW)- (3.6) 

We have the following CLT for functionals on interacting particle systems 
with finite range interactions, evaluated at a finite set of times. 

Theorem 3.2. (General CLT) Let d, M and H be as in Section 3.1. Suppose 
I C [0, 00). Suppose that for some 7 > 2, (3.3) holds for all t 6 I. Then for all 
s € I,t £ I, the sum 

a(s,t) := ^ Cov[H(CW),H(L z (^yj)} (3.7) 

converges absolutely, and for any B-valued sequence (A n ) n >i satisfying (3.4) 
and (3.5), as n — > 00 we have 

|A Il |- 1 Cov(5^(e- l/ ),^"(^"^)) - c(M) (3-8) 
and the finite- dimensional distributions of the process 

\K\- i/2 {sfr(it n n-Esfr(^ni tei (3.9) 

converge to those of a zero-mean Gaussian process with covariances given by 
(a(s,t),s,tE I). 
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When / = [0, r] for some r > 0, it is possible to extend the convergence in 
distribution of finite dimensional distributions to a functional CLT giving weak 
convergence in D[0, r] of the whole process (S^^f 71 '"), < t < r). Here D[0, r] 
is the space of right-continuous functions on [0, t] with left limits, equipped with 
the Skorohod topology (see [11]). One reason for an interest in such convergence 
is the study of weak convergence of hitting times (see [18]), although we do not 
pursue this further here. The functional CLT requires a stronger version of the 
moments condition (3.3), namely 

sup sup E[\H(tf ,v \tf)\ r ] < oo. (3.10) 

Aet30<t<r 

Theorem 3.3. Let d, J\f and H be as in Section 3.1. Let t > 0. Suppose 
(3.10) holds for some 7 > 8. Then with I = [0, r], the process defined by (3.9) 
converges weakly in D[0,t] to a zero-mean Gaussian process with covariance 
structure given by o~(s,t), s,t G [0, t]. Also, the limiting process admits a version 
with continuous sample paths. 

3.3. Remarks 

1. The limiting covariance in (3.8) does not depend on the choice of (A n ) n >i. 

2. We do not in general rule out the possibility that the limiting covariance 
in (3.8) might be zero, although in most applications this will not be the 
case. 

3. Theorem 3.2 can be extended to a CLT of measures, using results from 
[48]. Let 1Z denote the class of Ricmann measurable sets in R d (see [48] 
for definitions), fix Aq € 1Z with non-empty interior, and let A±, . . . ,Ak 
be Riemann measurable subsets of A$. Let (s n ,n > 1) be an increasing 
unbounded sequence of positive numbers, and let A i n := (s n Ai) P\Z d . Let 
ti,...,t)e be positive numbers. Then if (3.3) holds for some 7 > 2 and 
r > max(£i, . . . , t n ), we expect that 

{s- d/2 (s^atm - E[s^(ntr> v )\), 1 < i < fc ) 

converges, as n — > 00, to a multivariate centred normal with covariances 
given by a(ti, tj) x Lcb(|^ n Aj \) where Leb(-) denotes Lebesgue measure. 
It should be possible to prove this using theorem 2.1 of [48], but we have 
not written out the details. 

4. An alternative approach to the CLT would be via Stein's method, as used 
by Penrose and Sudbury ([49], theorem 7) for certain specific particle 
systems with finite local state space. It may be possible to adapt that 
approach to give a normal approximation in the more general setting of 
Section 2 (i.e., without assuming V is 1 d and with no stationarity as- 
sumptions), so long as we simplify H to a function H depending only on 
the local state (and not on neighbouring sites). That is, taking H to be 
a function from X to K, in the general setting of Section 2, one might 
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expect normal approximation for 

Again, we have not written out the details. 

5. These results are related to those given in Penrose [47] and in Doukhan 
et al. [18] for interacting particle systems with finite local state spaces. 
For other approaches to central limit theorems for interacting particle 
systems, sec Hollcy and Stroock [34]. See also De Masi and Presutti [16], 
Spohn [ 10], and Liggett ([40], page 39). These works are mainly concerned 
with cases where X is finite, and where there exists a unique invariant 
measure, which we do not assume here. 

6. The initial condition on the process namely independent and iden- 
tically distributed states at time zero, is perhaps more restrictive than 
one would like. We believe that it is possible to give similar results to 
Theorems 3.2 and 3.3 in the more general case of a spatially stationary 
initial distribution satisfying a strong mixing condition, if instead of the 
process £f we consider £ t , thereby making everything spatially stationary. 
Some other ways in which one might hope to extend the current work are 
described in Section 5.2. 

4. Spatial birth, death, migration and displacement processes 
4-1- General description 

Our setup is intended to be applicable to a large class of continuum particle 
systems. In these, the global state space S is the space of locally finite subsets 
of R d , d G N. 

Let R > 0, A > be fixed parameters; R denotes the range of interaction. 
Let Sr be the space of finite configurations of points in the ball Br centred at 
0. Given an element X of the state space S, for x £ R d let X\ x denote the set 
(— x + X) fl Br (an element of Sr). Consider the following possible transitions 
of the process when in state X. 

(i) Immigration ('birth'). Particles arrive as a space-time homogeneous Pois- 
son process of rate A. Incoming particles may be accepted or rejected. The 
probability of acceptance of an incoming particle at x is assumed to be a func- 
tion of the point process X\ x . If accepted, the incoming particle is placed at x 
and added to the configuration. 

(ii) Death. A particle at x disappears at rate S(X\ X ) where for each y G Sr, 
6(y) is a nonncgative number. 

(iii) Migration. A particle at x jumps at rate p(X\ x ), and when it jumps its 
new location x + y is distributed with distribution fi(X\ x ; dy), where for each 
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configuration y £ Sr, p(y) is a nonnegative number and /i(3^; •) is a probability 
measure on Br. 

To these ingredients, we can add the following embellishments: 

(a) Marked point process. Assume each particle carries a 'mark' in a measur- 
able space T (the mark space). Incoming particles have random marks, assumed 
independent and identically distributed according to some probability distribu- 
tion on T. The probability of acceptance of an particle at x may depend on its 
mark and those of the existing points in X\ x . A 'jump' (migration) of a particle 
may entail a change of the value of its mark as well as (or instead of) a change 
of location in R d . In the marked case, we modify S to be the set of locally finite 
subsets of M d x T; however, we shall often refer to subsets of M d x T more 
informally as 'marked subsets of R d '. 

(b) Synchronous updating. Whenever any one of the above types of event (im- 
migration, death, or migration) occurs, at a point x £ M. d , then this can cause 
some rearrangement, or change in the values of the marks, of the points within 
distance R of x (but without any change to the configuration at a distance more 
than R from x). It is assumed that this rearrangement occurs simultaneously 
with the immigration/migration/death event, and is determined by the nature 
of that event and the configuration of (marked) points within distance R of x. 
In the case of an immigration event at x, we allow the location and mark of 
the incoming particle itself to be affected by the existing points within distance 
R of x, so the particle may be placed at some point other than x but within 
distance R of x. 

Continuum processes of this type fit into the particle system picture of Sec- 
tions 2 and 3 as follows. Take the local state space X to consist of all finite 
subsets of [0, l) d x T (or more formally, finite point measures on [0, l) d x T). In 
the unmarked case we may think of T as having a single element. We partition 
R d into cubes C v := v + [0, l) d , v £ Z d , and the global state space S is identified 
with X z by identifying a locally finite marked point set X C M. d x T with an 
element rjx of X z obtained by setting 

T)x{v) = -v + (Xn(C v xT)), v£Z d , 

where for v £ Z d and y C R d x T we set v + y := {(v + x, 6) : (x, 9) £ y}. 
More informally, the local state of the system at site v £ Z d is the collection of 
marked points in C v . 

For d 6 Z d , we take Af v to be the set of w £ !* d such that the distance between 
the cubes C v and C v is at most R (and take N := Nq). When an immigration, 
death or migration event occurs at a point x £ C v , with or without synchronous 
updating, this may affect the states of other vertices, but only those in J\f v . 
Therefore, events at vertex v £ Z d affect only the states at vertices in M v , and 
this type of process fits into the setup of Section 2; in particular, (2.2) holds. 
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Moreover, rates are defined in a translation-invariant manner so that (3.1) holds. 

With the above identification of S with X zd , a Markov process in X 7 ^ corre- 
sponds to a Markov process in S. Given the local interactions (i.e., the param- 
eters of the immigration, death and migration rates), we can use Theorem 2.1 
to deduce existence of a Markov family processes (X t ,t > 0) in the state-space 
S with these local interactions. For bounded Borel A C M. d , let (<-f t A ,i > 0) 
denote the corresponding process when births, deaths and migrations taking 
place outside A are suppressed. 

A large class of stationary additive set functions on £ is given as follows. Let 
4> : T x Sji — > M be measurable and for bounded Borel Acl d , and X e S set 

Sf(X):= J2 <K9,(-x + X)n(B R xT)). 

The corresponding additive set function on X z is obtained by taking 
H^(vxW):= ]T (j>(j9,(-x + X)n(B R xT)), 

which depends on r\x only through T]x\jV- Given A £ B. let A denote the set 
L) ve AC v . We have the identity 

Sf(X) = S^(vx), 

so if satisfies the moment condition (3.3) for appropriate 7, we can use 
Theorems 3.1, 3.2 and 3.3 to obtain corresponding limit theorems for S^ n {X^ n ). 
The statements of results in Section 3 refer to processes 

0), taking- 
values in X Ijd , for finite A C 7L d . In the <S-valued process (xf' v ) corresponding 
to £^' v the immigration, death and migration events occur only in the bounded 
region A := U v ^aC v , but within this region the description of rates of these 
events is as before; we reserve the notation £ t for X z -valued processes and 
Xt for processes in S. Typically the physical system being modelled is in a 
bounded region, and Theorems 3.1, 3.2 and 3.3 provide information about the 
limit behaviour of the bounded-region process as this region becomes large. To 
apply these results, we need to specify v\ that is, we need to assume the initial 
configurations of (marked) points in each cube C z , z £ A are independent and 
identically distributed with some specified distribution v. For example, one could 
assume an initially empty configuration, i.e. assume v is concentrated on the 
configuration in Co with no points. 

To apply our general results, we need (2.4) to hold, i.e. we need a bound 
on the total jump rate at each site v <E Z d . With our description, the rate 
of immigration events per unit volume is bounded by A; however, if death or 
migration is allowed we need extra conditions to ensure that there is a uniform 
bound on the total rate at which these occur within a given cube C v . 

In some examples there is a hard-core constraint preventing particles from 
appearing at a distance less than e from each other, for some fixed e. In this 
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case there is a uniform bound on the number of particles that can coexist in C v , 
and provided the migration and death rates for each particle are bounded the 
total rate at which events happen in C v will be bounded. 

On the other hand, if the total number of potential particles coexisting in a 
cube C v is unbounded, for (2.4) to hold we need some stronger condition on the 
death and migration rates per particle than mere boundedness, since the jump 
rate at site v includes the sum of death and migration rates over all particles 
in C v . This could happen, for example, if the presence of a large number of 
particles nearby slows down the rate of updating of an individual particle in 
some suitable manner. 

Examples of the continuum particle system setup just described are numer- 
ous. Markov processes in S are discussed in [26; 27; 33; 52], and our results add 
to these works. Applications of these processes include the simulation of point 
processes in statistical literature (see e.g. [3; 23; 43]). In the case of deposition 
models, the 'state' of a particle in the set T could represent, for example, its 
shape, and/or its height above the surface. 

The following sections contain discussions of some examples. Our results 
(Theorems 2.1, 3.1 and 3.2) provide a systematic framework to approach ex- 
istence, spatial LLNs and spatial CLTs for these examples, although it should 
be pointed out that results along the lines of the LLN (Theorem 3.1) and our 
first CLT (Theorem 3.2, at least for onc-dimcnsional distributions) can also be 
achieved for some of these examples by ad-hoc modifications to the arguments 
in Penrose and Yukich [50]; indeed, this was shown by J.E. Yukich (personal 
communication 2001), who suggested many of these examples to the author. 

4-2. Random Sequential Adsorption and variants 

In the simplest version of Random Sequential Adsorption (RSA), particles do 
not carry marks and the only events are immigrations; there is no migration, 
death, or synchronous updating. An incoming particle is accepted with proba- 
bility 1 if no existing particle lies within unit distance of that particle, and 
otherwise. In effect, each particle lies at the centre of a hard sphere of unit di- 
ameter, and is accepted if its associated sphere does not overlap any previously 
accepted hard sphere. 

The model just described fits our framework, provided we take R > 1, and 
the existence result (Theorem 2.1) shows this process is well-defined in M. d . 
Our results give LLNs and CLTs for additive functionals of the form S£ on S, 
applied to the process (X^-), for which particles arrive only over the bounded 
region A C M. d . Note that in the basic model, qb is effectively a function on Sr 
rather than on T x Sr, because points are unmarked. Examples include the 
following: 

• If we take <f>i(X) = 1 for all X £ S R , then Sf^X) (for X £ S) simply 
counts the number of accepted points in A. 

• Let i?i be a parameter with 1 < R± < R, and let card (A 7 ) denote the 
number of elements of X. If (for X £ Sr) we take <j>2{X) = (l/2)card(A" fl 
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B^), then S£ 2 (X) (denned for X S S) counts the number of pairs of 
points within distance R\ of each other with both points in A, plus half 
the number of such pairs with one point in A. 

Generalizations of RSA include cooperative sequential adsorption (where the 
probability of acceptance is some function of the local configuration of points). 
RSA with the identical hard spheres of the basic model replaced by independent 
identically distributed randomly shaped hard objects (so particles carry marks, 
representing their shape), and RSA with desorption where particles leave (die) 
at a non-zero rate. RSA and variants are of widespread interest in applications, 
as discussed in [50], and references therein. 

Using Theorems 3.1 and 3.2, we may obtain LLNs and CLTs for S^ n (X t An ' v ) 
with i = 1 or i = 2. In the case of <pi, these have previously been given in 
[50]. Theorem 3.3 can be used to extend the CLT to a CLT in D[Q, oo) for 

(Sf;(X t A ^),t>0)- 

4-3. Continuum ballistic deposition 

Ballistic deposition (BD) models in the continuum have the following features. 
Particles are represented, typically, by unit diameter hard balls in d + 1 dimen- 
sions (usually d = 1 or d = 2). They fall vertically from above, sequentially 
at random, towards the adsorption surface R d (strictly speaking, R d x {0}). If 
a particle reaches the surface unhindered, it is attached to it. Details of what 
happens if a particle strikes another particle before reaching the surface depend 
on the model. 

In monolayer versions of BD, all accepted particles lie on the surface of the 
substrate, as in the basic RSA model, so necessarily some particles are rejected. 
In the version proposed for d = 1 by Solomon [59] and for d — 2 by Jullien and 
Meakin [37], on striking a previously deposited particle, the new particle rolls, 
following the path of steepest descent until it reaches a stable position. If this 
position touches the adsorption surface the particle is fixed there; otherwise the 
particle is rejected. It is assumed that the rolling occurs instantaneously. 

All accepted particles lie on the substrate, and so can be represented by points 
in M. d . The position of an accepted particle is a translate (or displacement) of 
the location in R d above which it originally comes in. There is a uniform bound 
i?2 on the size of the possible displacements (this is clear when d = 1, and 
a proof for d = 2 is given in [46]; see also [13]). Since the displacement and 
the decision on whether to accept an incoming particle at x are determined 
by the existing configuration within distance i?2 + 1 of x, this model fits into 
our framework. There are many other monolayer BD models satisfying this 
condition of uniformly bounded displacements. Again, the functionals S^(X) 
and S^ 2 (X) are of interest; for i?i = 1 the latter counts the number of touching 
pairs of particles in A. 

Senger et al. [56] describe many experimental results. Choi et al. [13] give 
both experimental and analytical results. Penrose [46] gives rigorous results on 
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the infinite input version of this model. 

In multilayer BD, a particle (ball) may attach itself to previously adsorbed 
particles instead of to the substrate. In the simplest form of continuum mul- 
tilayer BD, each particle falls vertically towards the substrate and as soon as 
it encounters either the substrate or another particle, it sticks (and remains in 
that place forever). 

In an alternative model of multilayer BD, when an incoming ball strikes a 
deposited ball, it does not stick but rolls, following the path of steepest descent, 
until it reaches a stable position. At this point it stops. 

In either of these versions of multilayer BD, each particle is accepted, and has 
a vertical displacement (or height) relative to the position it would occupy if it 
were to fall to the substrate unhindered by other particles. If we take T = [0, oo), 
this fits into our framework, and with the mark of a point representing the height 
above the substrate at which it is attached. In the case of multilayer BD with 
rolling, to apply our results we need a uniform bound on the horizontal distance 
a particle can roll. 

4-4- Deposition models with rearrangement 

The RSA and BD models can be made arguably more realistic if an incoming 
particle is allowed to cause some (instantaneous, finite range) displacement of 
existing particles. If we modify the RSA and monolayer BD models in this 
way, our general existence results, and limit theorems for functionals such as 
(<-f t A ) and S^ 2 (X^-) remain applicable, thereby adding to results in [50] which 
did not allow for synchronous updating, and to known experimental results for 
monolayer ballistic deposition with restructuring; see e.g. [37]. 

In the case of multilayer BD, such displacements can be viewed as 'avalanches'. 
Again, our general results remain applicable, provided we assume displacement 
has a finite range. Such an assumption is consistent with the assumption of high 
friction or bonding forces (see Onoda and Liniger [44]). 

Models of this type are related to the notion of random loose sphere packing 
(RLSP). This has long been of interest in Materials Science, since the classic 
experimental works of Bernal, Mason and Scott [9; 10; 55]. See Cumberland and 
Crawford [14] for a survey and over 300 references. While there is no one com- 
monly agreed upon mathematical model for RLSP, our multilayer BD model, 
with finite-range rolling and displacement allowed, has many of the features of 
the physical RLSP experiments. 

Among the additive set functions S£(X) of interest are those generated by 
the following choices of <j>: 

1. Total number of balls up to a specified height. Let i?3 > 0, and for X £ Sr 
let 4*3(0, X) be equal to 1 if 6 < R3, and otherwise. Then SjS (X) (defined 
for X G S) counts the number of points of X in A which have height at 
most i?3 above the substrate. If we take R3 = 1/2, then S^ 3 (X) gives the 
number of balls touching the substrate, since the balls have radius 1/2. 
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For large R3, S^ 3 (X) gives an estimate for the density of packed particles. 

Our results are applicable to give a LLN and CLT for S^ 3 n (Xf" n,v ) when 
(A n ,n > 1) satisfies (3.4) and (3.5). 

2. Total number of contacts ('coordination number'). Let 4>i(9, X) denote half 
the number of points (x, 6') E X such that ||(x, 9') - (0, 6)\\ = 1, where || • || 
denotes (d + l)-dimensional Euclidean distance. Then, modulo boundary 
effects, S^-(X) denotes the total number of contacts between particles in 
X. Our general results give us 

lim S^{X t A ^)/\A n \ = K{t) a.s. (4.1) 

with a corresponding CLT and functional CLT. 

3. Surface growth. In the multilayer BD process, possibly with rearrange- 
ment, the 'active zone' or 'interface' of the resulting agglomeration may 
be defined as the set of exposed particles, with a particle termed 'exposed' 
if there is a positive chance that the next incoming particle could strike 
it. Let $5(9, X) be the height 9 if a particle at (0,9) is exposed, and let 
<f>5(9,X) be zero otherwise. Then S£ (X ) is the total height of exposed 
particles over A. We can apply Theorems 3.1, 3.2 and 3.3 to get a LLN, 
CLT and functional CLT for S^(X t An ' u ). This adds to results in [51], 
where updating in the BD model is not assumed and the functional CLT 
is not considered. 

Estimates for coordination numbers are found in Bcrnal and Mason [10], 
Gervois et al [29] and Gotoh and Finney [30] (p. 202) and references 148-153 
of [14]. Jodrey and Tory [36] describe computer simulations of a loose sphere 
packing model. Non-rigorous work of Gervois et al [29] (p. 2124) suggests a 
coordination number of 6.1 ± 0.3 for a particular version of RLSP model. Gotoh 
and Finney [30] suggest a coordination number of 6.0. 

4-5. Off-Lattice interacting particles 

We indicate some continuum analogs of classic interacting particle systems on 
the lattice. 

a. ('voter model' I) T represents a set of two 'colours'. There is only immi- 
gration, with no death or migration. When a point arrives, it picks a point at 
random (within a distance R of it) and chooses its colour with probability p, 
< p < 1. If no existing point lies within distance R of the new point, its colour 
is selected at random. If 4> is the indicator of a particle having a given colour, 
then S^-(Xt) is the number of particles in A of that colour. 

b. ('Voter Model' II) When a point arrives, its adopts the colour of the nearest 
neighbour within distance R. If no such neighbour exists, it acquires its colour 
randomly. As the number of particles increases, the scale at which interactions 
takes place decreases. 

c. ('Exclusion process') Suppose particles do not arrive or depart, but jump 
in the state space. Suppose fl(dx) is a probability measure on Br, and A > 0, 
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e > arc constants. Suppose a particle at x £ M. d attempts to jump at rate A, 
and when an attempt is made to jump, the destination of the attempted jump 
lies in x + r with probability f3(T). The jump is successful if no previous particle 
lies within distance e of the destination. If the jump is unsuccessful the particle 
remains at x. 

d. ('Zero range process') A particle at x £ M. d jumps at rate A„ if there are n 
other particles distant at most e from x. When it jumps, its destination lies in 
x + T with probability (3(T). Here e and (3 are as in the preceding example. If the 
parameters A„ satisfy sup„(nA„) < oo then the bounded jump rate condition 
(2.4) holds. 

5. Further examples and discussion 
5.1. Discrete examples 

Not all applications of the general results of Sections 2 and 3 can be expressed 
in terms of marked point process in the continuum. Here we briefly discuss some 
cases where X is a discrete space. 

When X is finite, our limit theorems remain applicable. Indeed, if X is finite 
then any choice of H will satisfy (3.3) since H can take only finitely many values 
on a finite set. However, we concentrate here on examples where X is the infinite 



A model of multilayer lattice BD is be defined as follows. Given a vertex 
set V satisfying (2.2), assume particles arrive at each site v £ V as a Poisson 
process of intensity A. Each particle attaches itself to the system at a height 
which is 1 plus the maximum previous height of particles at sites in JV V . The 
particle system recording the maximum particle height at each site is a X v - 
valued process which jumps from state 77 to state rf at rate A for each vertex 
v £ V, where 



This type of model is of much interest in the physical sciences; see for example 
[5]. Mathematical studies [2; 57] have mostly been concerned with the one- 
dimensional case where V is Z or a sub- interval thereof. Our Theorem 2.1 shows 
the process is well-defined on more general graphs. Moreover, it can be shown 
that for V = Z d , when £0 = 0, and the additive set function is taken to 

be the sum ^2 v ga(£( v )) , for arbitrary k, then the moments condition (3.10) 
holds for all 7 and all r. 

A modification of the BD model is to allow surface relaxation (see [5]). Again, 
£t(w) denotes the height of the highest particle at site v at time t. But now an 
incoming particle at site v seeks the site w in J\f v with the lowest height £t (w) (in 
the event of a tie, select w at random from the minima) and locates itself there, 
adding 1 to the value of £t(w). Again, our results ensure existence of £t(w), and 
LLNs and CLTs for additive set functions of the form X™pa(£( w )) 



set {0, 1,2,.. .}. 
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Superficially related to these deposition models are 'sand-pile' type models, 
which have been much studied in recent years (see e.g. [17; 42]). However, they 
do not fit so easily into our framework, since there is typically no bound on the 
possible range of influence of a single incoming particle. 

5.2. Open problems 

Our general model has some limitations which restrict its scope; it would be of 
interest to remove some of these. 

One limitation already mentioned is the uniform bound on jump rates. This 
rules out direct application to a number of spatial immigration, death and mi- 
gration type models where there is no uniform bound on the total rates of 
migrations/deaths within a bounded region. 

Another limitation already mentioned is the restriction to finite-range inter- 
actions, which prevents our results being directly applicable to sand-pile type 
models, or to cooperative sequential adsorption models where the probability of 
acceptance of a particle depends on the whole existing configuration, not just 
within a finite range. 

In continuum spatial models, it is natural to try to allow for non-jumpwise 
motion. The simplest extension would be to allow for deterministic evolution of 
the system in between the jump-type events, with the deterministic evolution 
at one site affecting neighboring sites only through the jump rates. Such an 
extension could be sufficient to deal with models where the lifetime of a particle 
is not exponentially distributed; for a particular example of this sec the particle 
system in [58]. The theory of piecewise deterministic Markov processes [15; 
35] could be relevant here. In a more general extension, one might hope to 
allow piecewise deterministic motion of a particle. In this case, this motion is 
typically able to take the particle from one patch of space into a different patch 
of space (this feature does not arise in [58]), which is likely to make matters 
more complicated. 

6. Proof of Theorem 2.1 

For s £ F, and x £ X^ v , recall that a*(x, •) denotes the jump rate kernel of 
the Markov process on X^* with generator G*, and recall the definition (2.4) of 
Cmax as the supremum of all jump rates. As a first step, note that for x G X v , 
by (2.3) and (2.6) followed by (2.5) we have 

G v f(x) = [ ifZiv) - f v x {x\^))a* v {x\^,dy) 

= [ (f(x\v\y)-f(x))a* v (xW v ,dy). (6.1) 
Jx*v 

For finite A C V, let C(A) be the set of all bounded measurable functions 
/ : X v — ► M which depend only on co-ordinates in A, i.e. which satisfy 

f(x) = f(y), Vx,y e X v such that x(v) = y(v) Vv e A. 
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For finite A C V, let Ma := U veA M v . For / G C, let ||/|| := sup{/(x) : x G 
We shall use the following result to establish uniqueness of the Markovian 
family of processes with generator G. 

Lemma 6.1. Suppose A C V is finite and f G C(A). Then 

||G„/||<2c max ||/|| ) (6-2) 

G v f = 0, VveV\M A , (6.3) 

and 

G v f G C(AuM v ), VveMa (6.4) 

Proof. Since the total measure of a*(x|jv"„, •) is at most c max , (6.1) gives us 
(6.2). Also, if v G V\Ma and y G X^ v , then the restriction of x\v\y to A is the 
same as the restriction of x to A, so that f(x\v\y) = f(x), and hence G v f(x) = 
by (6.1). This yields (6.3). 

Now suppose v G Ma, and suppose x, x' G X v with x|ytujV„ = z'IaujV,,- Then 
a; agrees with x' on A4, so that the measures a*(x|v\/„, •) and a*(x'\j^ v , •) are 
identical. Moreover, x agrees with x' on ^4 so that fix 1 ) = f{x) and moreover 
for all y G we have that f(x\v\y) = f(x'\v\y). Hence, by (6.1) we have 

G v f(x) = [ (f(x'\v\y) - f(x'))a:(x'W v ,dy) = G v f(x'). 

Thus, G v f(x) depends on x only through x(w),w G A UM V , i.e. (6.4) holds. | 

On a suitable probability space, let (V v ,v G V) be a family of independent 
and identically distributed homogeneous Poisson point processes of intensity 
Cmax in [0, oo) x [0, 1]. Label the points of V v as ((Ti(v), Ui(v)), i = 1, 2, 3, . . .), 
with Ti(v) < T 2 (v) < T 3 (v) < ■ ■ ■ (with probability 1 the Ti(v) are distinct for 
each v). 

For A C V and v G V let M A denote the 2-ncighbourhood of A consisting of 
all vertices a graph distance at most 2 from some element of A in the graph V, 
and let M£ be the two-neighbourhood of v, i.e. let 

M+ := M Ma ; M+ := M+ } = M Mv ■ (6.5) 

By (2.2), the number of elements of M+ is at most 1 + D + D(D - 1) = 1 + D 2 . 
Define the point set P C V x [0, oo) by 

P := U„ ey ({(«, 0)} U {(v, Ti(v)) : i > 1}) . (6.6) 

Make into the vertex set of an (infinite) oriented graph Q by putting in an 
oriented edge (u, T) -> (v, T') whenever v G M£ and T < T'. 

For w, v G V, let us say that w affects v before time t if there exists a 
(directed) path in the oriented graph that starts at (w, 0) and ends at some 
Poisson point (v,T) with T <t. Let E t (w,v) denote the event that w affects v 
before time t. Let dist(-, •) denote graph distance. 
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Lemma 6.2. There exists a constant 5 > 0, such that for all v G V, and all 

n G N, 



u 



[E Sn (w,v)} 



w £ V : dist ( w ,v ) > 2n 



u 



[^(^w)] 



u; £ V : d is t ( w , v ) > 2 n 



< 2~ 



< 2" 



(6.7) 



(6.8) 



Proof. Let be exponentially distributed with mean l/c max , and let S n be 
the sum of n independent copies of W. Let 9 > 0. Then by similar arguments 
to those used in the proof of Lemma 5.1 of [47], the probability on the left side 
of (6.7) is bounded by {D 2 ) n P[S n < Sn], and hence by D 2n e e5n (E[ e - ew ]) n . 

Choose 9 so that E[e~ ew ] < 1/(4D 2 ), and then choose d > so that e es < 2. 
Then the preceding bound is at most (1/2)™, completing the proof of (6.7). The 
proof of (6.8) is similar. | 

Corollary 6.1. Let v G V and t > 0. Then (a) wit/i probability 1, v affects only 
finitely many w G V before time t, and (b) with probability 1, only finitely many 
w G V affect v before time t. 

Proof. Let F n be the event that there exists w € V with dist(u;, v) > 2n such 
that v affects w before time t. Then P[F n ] — > as n — > oo by (6.8), and this 
implies part (a). The proof of (b) using (6.7) is similar. | 



We partition the vertex set V of G into 'generations' by putting each point 
(v, 0),v G V in generation and saying a point (v,T), T G (0,oo) of T 3 is in 
the fcth generation if the longest directed path from any point (w, 0) to (v, T) is 
of length k (with length measured by graph distance). By Corollary 6.1, this is 
indeed a partition with probability 1. 

Recall that a*(x) denotes the total measure of the jump rate kernel, and 
that a* (a;) < c max . Let S x (-) denote the unit point mass at x and define the 
probability kernel 

^ {xr)= <M + {^-at{x))8 x {.) (g9) 

Cmax 

For each v G V, there is a measurable function ip v : x [0, 1] i— > X^ u 

with the property that if d is uniformly distributed over [0, 1], then ij) v (x, z?) has 
distribution fx v (x, •) for each x G X^ v . The existence of such a ^ follows by 
Kallcnbcrg ([38], page 56, lemma 3.22), since [i v (x, B) is a probability kernel on 
X*«, which is a Polish space, and hence a Borel space (see [38] theorem A1.2, 
or [24], p. 419, Proposition 20). 

For v £ V, set Tq(v) = and list the arrival times Tj(u), u G Af v in increasing 
order as T*(w), T 2 *(u), T 3 *(w), .... 

Given x G X v , construct the process (£f,i > 0) as follows. First, for all 
uGV assume that £f (u) is a right-continuous function of £ that is constant in 
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between arrival times T*(v),i > 1. In other words, for all v G V and t > 0, if 
< * < T*(v) then set £f(v) = £?.. («). Set = for all sel/. 

Suppose inductively that £f.(i>) is defined for each (v,T) in generation k of P. 

Suppose arrival (v,Tj(v)) is in generation fc + 1. Set T = Tj(v). Dehne the 
new state of the process in M v at time T by 

@k=iM@_k,.tfj(«))- ( 6 - 10 ) 

This is well-defined because for each w G there exists i = i(w,T) such that 
T = T*(w), and (w,T*_ 1 (w)) is in generation ft at the latest, so that &(w) 
is already defined for t in the interval [T*_ 1 (w),T) and hence £j,_(w) is well- 
defined. 

Lemma 6.3. There exists a filtration (Qt,t > 0) with respect to which the 
family of processes ((^,t >0),xG X v ) described above is a Markovian family 
of processes in state-space X . The associated family of transition distributions 
induces a semigroup of operators (Pt,t > 0) with generator G given on C by 
(2.7); that is, with this Pt and G, (2.1) holds with uniform convergence for all 
feC. 

Proof. Let Gt denote the cr-field generated by all Poisson arrivals before time 
t at all sites, i.e. let Qt denote the smallest cr-field with respect to which each of 
the random variables of the form V V (B), with v £ V and B a Borel subset of 
[0,i] x [0, 1], is measurable. Here V V (A) denotes the number of points in A of 
the Poisson process V v ■ 

By the construction and standard properties of the Poisson process, the pro- 
cesses (£ x ,x € X v ) form a Markov family of processes, adapted to (Qt,t > 0). 
Sample paths are right-continuous with left limits because we use the product 
topology and there are no accumulations of Poisson arrivals at any site v G V . 

Let (Pt,t > 0) be the semigroup associated with the Markov family (£?,x G 
X v ). Let / G C and let x G X v . Then 

P t f(x)-f(x)=E[f(Q)-f(x)\. 

Let x G X v . Choose finite A C V such that / G C(A). Then /(£f) = f(x) 
unless min{Ti(w),w G Ma} < t. 

Let be the event that that T\(v) < t for two or more v G N^, and let 
F^ 2) be the event that T 2 (v) < t for some v G Ma- Let F t := F t {1) U F^ 2) . Also, 
for v G Ma, let F$ be the event that T x (v) < t < T 2 (v) and T x (w) > t for 
w£M£\ {v}. Then 

Eif(&) - f(x)\F$] = [ (f(x\v\y) - f(x))f, v (xW v ,dy) 
= c~L / (f(x\v\y) - f(x))a* v (x\ Nv ,dy), 

so that by (6.1), 

E[f(m-m\F$\ = c-i x G v f(x). 
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Hence, 

Em) m 

= P[F t ]E[m)-f(x)\F t ] + Yl P[F^}E[f(^)-f(x)\F t {3) ] 

= P[F t ]E[f(Q) - f(x)\F t ] + P[F$]^G v f(x). (6.11) 

v 6 A/a 

Using (6.11), then (6.3) and then (6.2), we obtain 

|t- 1 ( J B/(^)-/( 2 ;))-G/(.T)| 

= r i P[F t ]E[f(e t )-f{x)\F t ]+ J2 <r l <£*p[F$!] - 



v£Af A 



< 2||/|| ( r l P[F t ] + (t-'c-^PiF^] - l) j . 



(6.12) 



As t i we have that P[F t ] = o(t) and P[F$] = c max t + o(t): moreover P[F$] 
is the same for all v <E Ma- Hence the bound in (6.12) is independent of x and 
tends to as t J, 0. Thus the generator of the semigroup {Pt,t > 0) is G on 
functions /eC.| 

The next step is to show that the transition semigroup (Pt,t > 0) satisfies 
(2.8). For finite A C V, to get a process with generator G A rather than G, we 
need to 'switch off' all jumps taking places at sites outside A. To do this, set 

V A :=Tn(Ax [0,oo)), 

and partition V a into generations in precisely the same manner as for V ■ Then 
using the updating rule (6.10) but now restricting attention to arrival times 
Tj(y),v € A, we obtain a pure jump process (£ t ' x f t > 0) with generator Ga- 
Note that we use the same family of Poisson processes {V v , v £ V ) in construct- 
ing (£ t ,a: , t > 0) and (^f , t > 0), regardless of A. 

For each («, t) € V x [0, oo), let the 'cluster' C^ jt be the set oiw £V such that 
w affects z before time t for some z G A/J" (recall A/"+ is the 2-neighbourhood 
of v). By Corollary 6.1 this cluster is almost surely finite. Its significance lies 
in the following result, which will be used again in proving limit theorems later 
on. 

Lemma 6.4. Suppose (v, t) G V x [0, oo) 7 and x £ X v , and A C V is finite. If 
C v , t C A, then £?' x {w) = Q{w) for all w E M v . 

Proof. All differences between (£,f' x ,s > 0) and s > 0) are due to Poisson 
arrivals in (V \ A) x [0, oo) x [0, 1] which contribute to changes in (£f , s > 0) 
but not to changes in (^' x ,s > 0). However, by the very definition of C V) t, 
since C v t C A the influence of these Poisson arrivals does not propagate to any 
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vertex in by time t, and therefore £f |jv„ = £,t' X \N v - 



Lemma 6.5. The transition semigroup (Pt,t > 0) of the Markovian family of 
processes ((^f , £ > 0),x £ X v ) defined above satisfies (2.8). 

Proof. Suppose / £ C and choose finite A C V such that / £ C(A). Suppose 
iW m , m > 1) is a sequence of finite subsets of V satisfying liminfm^oo W m = V. 
Then for v £ A, by Lemma 6.4 we have £ t m,x {v) = £f(v) for large enough m, 
almost surely. Hence, /(£ t m,x ) = f(Q) for large enough m, almost surely. Hence 
by the dominated convergence theorem, 



Ptf(x)=E[f(g)]= lim E[f{aY m ' X )] 



lim P^ m (x) 



Proof of Theorem 2.1. In view of Lemmas 6.3 and 6.5, it remains only to 
show that the transition semigroup (Pt,t > 0) is the only one with generator 
G on functions in C. We do this directly, without appealing to generalities on 
resolvents, in part because the state space is not locally compact in general. 
Let f £ C, and choose a finite set A C V such that / £ C(A). Write Af£ for 
Uugx-A/^- By Lemma 6.1, 



Gf=J2 G -f e c - 



(6.13) 



For t > and ft. > 0, we have 



Ih-^Pt+uf - P t f) - P t Gf\\ = 



P 



Phf - f 



-Gf 



as h i 0, 



because by assumption G is the generator, defined on / £ C by assumption so 
that (2.1) holds with uniform convergence, and Pt is a contraction. Similarly, 
for t > and < h < t, 



\h~\P t f - P t - h f) - P t Gf\ 



Pt-h 



Phf - f 



Gf 



as h | 0. 



Hence for any x £ X v , we have f t P t f{x) = P t Gf{x). Since Gf £ C by (6.13), 
for integer n > we can repeat the preceding argument n times to deduce that 



Ptf(x) = P t G n f{x) 



Hence by Taylor's theorem with Lagrange remainder, 



Ptf{x) 



G j f(x)t j ) P u G n+1 f(x)t n+1 



some u £ (0,t). (6.14) 
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By twice applying Lemma 6.1, we have that for v e Ma and w 6 V, if 
w NauN v then G w G v f = 0, whilst in any case G w G v f <G C(AUJV V U W„,). 
Hence, if w N~a\j{v} ( as defined in (6.5)) then G w G v f = 0, whilst in any 
case G w G v f € C(AL)Aft VtW \). Repeating the same argument, an induction on n 
shows that G Vn G Vn _ 1 ■ ■ ■ G Vl f = except when each vt lies in ■Af^u{v 1 «• i}" 
Hence, iterating (6.13) yields 

G n f(x)= J2 E E G Vn G Vn ^---G v J(x). (6.15) 

Recall the definition (2.2) of -D, and let | ■ | denote cardinality. The number 
of n-tuples (v\ . . . , v n ) £ V n such that each Vi lies in M^ u ^ Vl v . ^ is bounded 
by the expression 

\A\(\A\ + 1)(\A\ + 2) • • • (\A\ + n- 1)(1 + D 2 ) n 
< 2™ max(|A|, 1) max(|A|, 2) • • • max(|A|, n){\ + D 2 ) n 

< n\\A\^ A h n {l + D 2 ) n . 

Moreover, by (6.2), for all Vi,v 2 , ...,v n we have \\G Vn ■ ■ ■ G V J\\ < (2c max )"||/||. 
Thus, setting Kq := 4c max (l+Z) 2 ), by (6.15) and the fact that P u is a contraction 
we have that 

\\P u G n f\\ < \\G n f\\ < n!|^|l A l/f "||/||, u > 0. 
Thus for t < 1/Kq, the remainder term in (6.14) tends to zero as n — > oo, so 

oo 

Ptf{x) - J2 > J(;J !'■>■■ t < V^o. 

Thus, for / 6 C and 4 < 1/Kq, we have an explicit expression for Pt(f) in terms 
of G. Hence, if (A)t>o is another semigroup, associated with a Markov family 
of processes in X v , with generator G on C, we must have Pt(f) = Pt(f) for 
all t < 1/Kq. By the Monotone Class Theorem [62], this identity extends to all 
bounded measurable functions on X v , and by the semigroup property it then 
extends to all t > 0. This completes the proof of uniqueness. | 

7. Proof of general limit theorems 

To prove Theorems 3.1 and 3.2 we shall apply results from [47]. Recall that 
V = Z d for these results and B denotes the collection of finite subsets of Z d . 
Given a probability distribution v on I, let the process {£," ,t > 0) and (for 
A G B) the process (£ { ' v ,t > 0) be as described in Section 3.1. That is, 
(respectively £ A '") has generator G (respectively G A ) with initial distribution 
given by the product measure with all marginals equal to v. 
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We assume for our proofs that the process , t > 0) is constructed in terms 
of the family (V v ,v G Z d ) of Poisson processes in [0, oo) x (0, 1), as described in 
Section 6, in particular the description leading up to (6.10), with initial values 
(£o v ^ gi yen by a family of independent X-valued variables v € Z d ) 
with common distribution independent of the Poisson processes V v ,v G Z d . 

Likewise we assume the process {^f'^jt > 0) is constructed using the same 
family of Poisson processes, as described just before the statement of Lemma 
6.4, and using the same initial values (Cv,v G 7L d \ Thus we couple the processes 

and £. so the initial values are identical £q = £ ' v and the same Poisson 
processes (V v ,v G Z d ) are used in both cases (but in the case of £^ the effects 
of the Poisson processes V v for v ^ A are 'switched off'). 

Recall from Section 6 that for (v, t) G Z d x [0, oo) the 'cluster' C V: t is the set 
of u G Z d such that u affects v before time t for some v G AfJ ', and is almost 
surely finite. Since with our coupling wc have £q = £jf the following result is 
immediate from Lemma 6.4. 

Lemma 7.1. Suppose (v,t) G Z d x [0, oo) and suppose 
Cu,t C A, then £f' u (w) = £^ (w) for all w G JV V . 

Let i > 0. Given A e B and n e 4, i > 0, define F„(A, i) by 



if ue A\A° 



If (Wa) = A, we have 



ue(A0i) ueA^i 

In the terminology of Section 3 of [47], for any t > the random variables 
Y z (A,t), defined for z G A, A G B, form a stationary B-indexed summand with 
respect to the i.i.d. family ((V z , Cz), z G Z d ). To prove Theorems 3.1 and 3.2, we 
shall use Theorem 3.1 of [47]. We need to check the conditions of this result. 

Proof of Theorem 3.1. If (A n ,n > 1) is a B- valued sequence with liminf A„ = 
Z d , then with probability 1, A n includes the cluster Co tT for all large enough 
n, and hence by Lemma 7.1, ^ n ' v {v) = £"(v) for all but finitely many values 
of n, and all v 6 Af. Hence, Y {A n ,r) tends to Yb, T := H(£_"\_\f) as n — > oo, 
almost surely. Moreover, the variables Yo(A n , r) are uniformly integrable by the 
assumption that (3.3) holds for some 7 > 1, and therefore the convergence of 
Yo(A n ,T) to Yo.t extends to convergence in L 1 . 

Thus, all the conditions for the first part of Theorem 3.1 of [47] hold (note 
that \A n \/\Af An \ -> 1 by (3.4)), and by (3.3) of [47] we obtain (3.6). | 

We now work towards a proof of Theorem 3.2. Given t > 0, v G Z d , and 
A G B, define the random variables 

Y Vit := H(L V ($W V )); Y$ := H(L v (^W v )). (7.2) 
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We now give some bounds on covariances. For v E Z d and A C Z d with A ^ 0, 
set dist(i>, A) := min{\v — w\ : w E ^4}. Also, recall from Section 3.1 that B r is 
the closed Euclidean ball centred at the origin in M. d , and for r > and v E Z , 
write B r (y) for v + B r . 

Lemma 7.2. Let s > 0,i > and assume the moments condition (3.3) holds 
for t E {s,t}, for some 7 > 2. Then there is a constant K such that for all 
v E Z d , w E 1 d and A E B we have 

\Cov(Y VtS ,Y Wtt )\ <Kexp(-\v-w\/K) (7.3) 

and this also holds ifY VfS is replaced by Yy^ and/or Y Wyt is replaced by Y^^, 
and also 

\Cov{Y v W - Y Vtt ,Y$)\ < Xexp(-dist( w , A c )) (7.4) 

and this still holds ifY^ A J is replaced by Y w t . 

Proof. Let v <EZ d ,w E Z d . Let r := (\w-v\/3). Let Y v>s := 3^, 8 l{c„,,CB r (t;)}- 
Then Y VtS is determined by ((P Z ,Q), z 6 J Br(w)nZ d ), so that Gov(Y VyS ,Y Wyt ) = 0. 
Choose 7 > 2 such that (3.3) holds for r E By Holder's inequality and 

Lemma 6.2, there exists a constant K (independent of v,w) such that 

E[(Y ViS - Y VyS ) 2 } = E[Y 2 s l {Cv y3 \B r (v)^9}] 

< (E[Y^ S ]) 2 ^P[C V , S \ B r (v) ? %]^)h 

< Kcxp(—r/K) 

and similarly, E[(Y w , t — Y UKt ) 2 ] < Kexp(—r/K) so that using the Cauchy- 
Schwarz inequality and the moments condition (3.3), we obtain for another 
constant (also denoted K) that 

\Cov(Y ViS ,Y W!t )\ < \Cov(Y ViS -Y ViS ,Y W!t )\ + \Coy(Y V!S ,Y W!t -Y W!t )\ 
< Kexp(-r/K) 

so that (7.3) holds. The proof is identical with Y VfS replaced by Y^ or Y w _t 
replaced by Y^J . 

Turning to (7.4), put r' = dist(v,A c ) - 1. Then Y v>3 = Y$ unless C v , t \ 
B r /(v) 0, so by Holder's inequality, (3.3) and Lemma 6.2, 

E[(Y$ - Y VyS ) 2 ] < E[(Y V W - Y VtS y] 2 ^P[C v , t \ B r ,{v) ± ^- 2 )h 
< Kexp(-r'/K), 

and combining this with the Cauchy-Schwarz inequality and using (3.3), we ob- 
tain (7.4). | 

We now give a limit for covariances, which is part of the statement of Theorem 
3.2. 
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Lemma 7.3. Let s > and t > 0. Suppose for some 7 > 2 that (3.3) holds for 
t G {s,t}. Suppose (^4 ra )n>i is a B-valued sequence satisfying (3.4) and (3.5). 
Define o~(s,t) by (3.7). Then asn-> 00, 



(7.5) 



Proof. Observe that ((Y VtS ,Y v j),v G Z d ) (defined by (7.2)) form a stationary 
random field. For r > 0, let A r n be the set of sites v G 1 d such that B r (v) flZ d C 
A„ . Then for any constant r > 0, 

\A n \- 1 Cov(Sfr(£),Sfr(%)) = \A n \- 1 E Cov(Y v , s ,Y w>t ) 

veA n w£A n 



E E cov^,-.^) 



(7.6) 



Every v G A n \ A r n lies in B r (v') for some 1/ G d ex t(A n ), so that |^4„ \ A r n \ < 
\B r n Z rf | x \d exk (A n )\, and so by (3.4), as n -> 00 we have |A„ \ A£|/|A„| -> 
so that |A^|/|A„| — > 1. Also, by (7.3), for some constant K\, 



\A ri 



J2 E CorWv^Y^t) 

>eA n \A^ w£A„ 



^tfilA^IAA 4ll-o. (7.7) 



Also, 



limsupl^r 1 E E |Cov(^, s ,F w>t )|< ^ |GOT(y ,.,r a>t )|:=^ 1 (r), 

w£Z d \A n zGZ d \B r 

where h\(r) — > as r — > 00 by (7.3). Combining this with (7.7) in (7.6), we 
obtain the limit 



lim \A n \- 1 Cov(Sfr(H v s ),Sfr(%))=a(s,t). 



(7.8) 



Next, writing Y^" for Yj t , we approximate the non-stationary random field 

((Y$,Y$),v G Z d ) by the stationary random field ((Y ViS ,Y Vtt ),v G 1 d ). We 
write 

Cov(s£" (^),S^ (^)) ~ Cov(Sfr (£), ^" (#)) 

= E E (Cov(yW,yW) - Cov(y;, s ,r^ t )) 



f E E OovcyW-^yW)] 
+ E E cov^.y^-y^t). 



(7.9) 



vGA n weA„ 
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By Lemma 7.2, there is a constant K 2 such that 

\Cav(Y$ - Y V ,„Y$)\ < K 2 exp(- max(dist(«, A c n ), \v - W |)/A' 2 ) (7.10) 
and from this we may deduce that for v £ A r n we have 
\Cov(X$ -Y Vt .,Y$)\ 

< J2 K 2 cM-\z\/K 2 )+ J2 K 2 exp{-r/K 2 ) 
=■■ h 2 (r) 

where h 2 (r) — ► as r — > 00. Also, for any r we have by (7.10) that for some 
constant A3, 

i^r 1 E E i Cov ( y ^ - Y «,» Y 3)\ < K 3 \A n nA n \ <i - . 

Thus, by splitting the v sum into v £ A r n and u £ A„ \ A r n we obtain that 
hmsupKr 1 ]T £ |Cov(yW-n, S) y^)|<ft 2 (r) 

and since h 2 (r) — > as r — > 00, this shows that 

hm E E \Cov(Y$-Y v ,.,Y<$)\=0. 

v£A„ w£A n 

A similar argument shows that 

lim M^l E E |Cov(^ s ,y^-^ t )| = o. 

Combining the last two limiting expressions, we obtain from (7.9) that 

and combining this with (7.8) we obtain (7.5). | 

Proof of Theorem 3.2. As usual, for any random variable £- and any p > 1 
we write := A[|^| p ] 1 / p when the expectation exists. Recall the definition 
of Y Z (A, s) at (7.1). Define Y z (A,t) similarly but with (Vo,(o) rcsampled (i.e., 
replaced by an independent copy of (Vo, Co))- 

For t > 0, let Dt denote the set of v £ Z d such that there exists w £ Af v such 
that affects w by time t. If we resample (To, Co) (leaving ((V v , (v), v £ Z d \{0}) 
unchanged), then the set D t is unchanged; this is because because (0,0) is 
included in V by (6.6) and any point (v, T) £ V with v ^ that can be reached 
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by a path starting at (0, 0) can be reached by such a path that does not use 
any point (0,T) with T > (see Section 6). By part (a) of Corollary 6.1, the 
set D t is almost surely finite. Also, for any A £ B, if v ^ D t then we have 
Y V (A, s) = Y£(A, s) for all s<t. 

We assume there exists 7 > 2 such that (3.3) holds for all t £ I. Choose 
7' £ (2,7). For any A € B we have by the Minkowski and Holder inequalities 
that 



veD t nA 



Y v (A,t) 



1' 



J2y v (A, t)l {veD 



} < \\ Y v(A,t)i {veDt} \w 

< J2 \\Y v (A,t)\\ 7 P[v £ £> t ](V7')-(i/7) 

vez d 

and by (3.3) along with Lemma 6.2, this is bounded by a constant that does 
not depend on A. 

Let k £ N. Let (t\, t^, ■ ■ ■ , tk) £ / with distinct components, and let (61, ... , &&) G 
M fe . Consider the linear combination R n := J2i=i hSff n (^ t .""' / ). Using (3.3) for 
7 > 2, we can apply Theorem 3.1 of [47] to deduce that there exists a 2 > 0, 
dependent on (&j,ii)f =1 , such that as n — > 00, 

LA^Var^) o- 2 , l^l" 1 / 2 ^-^) ^iV(0,a 2 ). (7.11) 

By the first part of (7.11), and Lemma 7.3, 

k k 

i=i j=i 

with er(s,t) given by (3.7). Set 

S^'") := S^'") - BS^'")- (7.12) 

By the second part of (7.11) and the Cramer- Wold device, the finite-dimensional 
distributions of the process {\A n \~ 1 / 2 S^ 1 n (^f' n,l/ ),t £ I) converge to those of a 
Gaussian process with covariance function (cr(s, t),s,t £ I), and this completes 
the proof. | 

Proof of Theorem 3.3. Again we use the notation S^^f'") given by (7.12), 
for the centred version of S$(£t )• We have shown convergence of the finite- 
dimensional distributions of the sequence of processes (\A n \~ x t 2 S H n (£,. n ' v ), n > 
1), and it remains to prove tightness in £)[0, t] of this sequence of processes. 
We follow the standard procedure of demonstrating tightness by estimating 
moments of increments, as in Billingsley [11], Theorem 15.6. 

Let < s < t < t with t — s < 1/2. Changing notation from the previous 
proof, for v £ 1i d we set 



Y v (s,t) 



H(L V (^WJ) - H(L V (^>»W V )). 
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Then 

(Y v (s,t)-EY v (s,t))\ 

4 

= E E E E El[(Y Vi (s,t)-EY Vi (s,t)), (7.13) 

vi£A n v 2 eA n v 3 eA„ v±£A„ i=l 

For v := (ui, . . . , Vk) E (Z d ) k , each (1 < i < k) is an element of Z d which we 
shall refer to as a d- component of v. We shall write Vi € v as shorthand for l Vi 
is a d-component of v'. 

At this point our argument is loosely inspired by the proof of Proposition 
8.7 of [18]. Given v = (vi, i>2, ■ • • ,v q ) E (Z d ) 9 and r > 0, let us say that v is 
r-connected if the graph with vertex set {1, . . . , q} and an edge between each 
pair such that \\vi — Wj||oo < r is connected. Define the gap of v to be the 
smallest r such that v is r-connected. Let G r {q,n) be the set of v S {A n ) q with 
gap r. 

For v := (vi, v 2 , v 3 , v A ) E (Z d ) 4 , and y 6 (Z d ) k and z E (Z d ) 4 ~ fe with k E 
{1, 2, 3}, we say y and z form a splitting of v if there exists nonempty proper 
subset J of {1, 2, 3, 4} such that y = i E J) and z = (u,, i € {1, 2, 3, 4} \ J). 
Define the distance between y and z to be min{||y — z||oo} ; where the minimum 
is over all c?-componcnts y E y and z E z. 

Suppose v = . . . , U4) £ G r (4, n). Then we assert that v admits a splitting 
(y, z) such that (i) the distance between y and and z equals r, and (ii) both y 
and z are r-connected. To see this, let the edges {i, j} of the complete graph on 
vertex set {1, 2, 3, 4} be assigned weights Wij = \\vi — vj H^, and let these edges 
be removed one by one, in order of decreasing weight, until the graph becomes 
disconnected. Then the last removed edge has weight r, and after removing this 
edge the graph has precisely two components which determine a splitting with 
properties (i) and (ii) above. 

Given v = (v\, . . . ,v 4 ) E G r (4,n), let (y(v),z(v)) be a splitting of v with 
properties (i) and (ii) of the previous paragraph. Let the choice of (y(v),z(v)), 
out of all such splittings, be made by an arbitrary deterministic rule. For y E 
(R d ) k set 



E 



A n ,v< 



Hy ~l[(Y y (s,t)-EY y (s,t)). 



(7.14) 
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Then 

4 oo 

E E E E Ei[(Y Vt ( S ,t)-EY Vt ( S ,t)) = J2 E £ i n vi 

v 2 GA n v 3 £A n V4€A n i=l r=0 vGG T -(4,n) 

00 

= E E E[IL y(v) U x(v) ] 

r-0 V £G r (4,n) 

D \ OO 

J2 E Cov(n y(v) ,n z(v) ) + £ ^[n y( v)]^[n z(v) ]. (7.15) 

^r=0 v£G r (4,n) / r=OveG r (4,n) 

If v G G r (4, n), then y(v) and z(v) are both r-connected. Hence, with | ■ | 
denoting cardinality, there exists a constant K3 such that 

\G r (k,n)\ < K^A^r*- 1 ^)*- 2 = Kz\A n \r d ( k -V-\ (7.16) 

Next, we assert that there exists a constant K4 such that for v £ G r (k,n), 

|Cov(n y(v) ,n z(v) )| < K 4 cxp(-r/K 4 ). (7.17) 

The proof of this assertion is similar to that of (7.3). Put 



: = II ((Y y ( S ,t) - EY y ( S ,t))l {Cy t CB r 



,/ 3 (y)} 

and choose 7 > 6 such that (3.10) holds. Then by Holder's inequality and 
Lemma 6.2, 

£[(H y(v) - H y(v) ) 2 ] < (E[\U y{v) \^P[U yey{v) {C ytt \ (y + B r/3 ) ± %}]^h 
< Kex-p(-r/K) (7.18) 

and likewise for z(v). Since n y ( v ) and n z ( v ) arc independent, we obtain 

|Cov(n y(v) ,n z(v) )| < |Cov(n y(v) -n y(v) ,n z(v) )| + |Cov(n y(v) ,n z(v) -n z(v) ))| 

and using (7.18), along with the analogous expression for z(v), and (3.10) and 
Holder's inequality gives us (7.17) as asserted. 

Next we estimate n y ( v ) in terms of (i — s) for small t — s. Recall that v = 
(i>i, . . . , W4) and y(v) = (vt, i G J) for some nonempty proper subset J = J(v) 
of {1,2,3,4}. Expanding the product in (7.14), we obtain 

n y(v)= E f II (-E[Y Vl (s,t)})) x l[Y Vj ( S: t), (7.19) 

where the product over the empty set is 1. Denote by II ji the last factor in 
(7.19), i.e. n,/' := Ylj £ jiY v .(s,t). Then II j< = unless at least one of the 
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Poisson processes V w ,w G Ui e j>Af^ has an arrival between time s and time 
t. Hence, if k denotes the number of elements of J then if J' ^ there are 
constants K$,K' 5 such that 

E[tl 2 j,} < E[\Uj,\i' k f k /iP[tij, ± O^-W-r) < K 5 (t - sf-^h) 

so that ElP (v) < K' 5 (t-s) l -V k l~<\ and similarly E^ 2 ,,] < K'^t-s) 1 -^ 4 -^/^, 
so that by the Cauchy-Schwarz inequality there is a constant K§ such that 



Cov(n y(v) ,n z(v) ) <Kg(t-s) 



l-(4/ 7 ) 



(7.20) 



Choose a > 1/2 with a < 1 — (4/7) (this is possible because we assume 7 > 8). 
By (7.16), (7.17) and (7.20), there are constants Kq and K7 such that the first 
term in the last line of (7.15) satisfies 



E E Cov(n y(v) ,n z(v) ) 

r=0 veG r (4,n) 



< K 6 \A n \ 



r 3d 1 min(e 



-r/if 6 



r>0 



53 r^-^t-s) 1 "^) 

, 0<r<(l-(4/ 7 ))if 6 | log(t-s) 



E 

r>(l-(4/ 7 ))JC 6 |log(t-«)| 

< K 7 \A n \(t-s) a . 



r 3d-l e -r/K B 



(7.21) 



To estimate the second term in the last line of (7.15), note first that by the same 
argument as for (7.21), there is a constant Kg, such that 

00 

E E \m y }\<K S \A n \(t- s r. 

r=0 y£G r (2,n) 

Then observe that i?[n y ( v )]i?[n z ( v )] = unless both y(v) and z(v) have two 
d-componcnts, and hence 



53 53 \E[IL Hv) }E[IL <v) }\< 53 53 \E[U y 

r=0v£G r (4,n) \r=0 y<EG r (2,ri) 



(7.22) 



Combining (7.21) and (7.22) in (7.15) we find that for some constant denoted 
Kq, the right hand side of (7.13) is bounded by 



K 9 (\A n \(t - s) a + \A n \ 2 (t - s) 2a ), 
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and hence by (7.13), 



E 



(7.23) 



We now set U n {t) := \A n \- l / 2 S^{^ n ^). Given 5 > with 8 < 1/2, by (7.23) 
and the Cauchy-Schwarz inequality we can find no = no (5) such that for n > no 
and 0<t 1 <t<t 2 <T with i 2 — h < 8, 

E[\U n (t) - U n (h)\ 2 \U n (t 2 ) - U n (t)\ 2 } < 2K 9 S 2a =: (F(8)) 2a . 

We slightly modify the argument of proof of Theorem 15.6 of Billingsley [11], 
using notation wf from (14.2) of [11] and w"(U n ,6) from (14.44) of [11]. Given 
e > and 8 > 0, choose no (6) as above. By following pages 129-130 of [11] we 
obtain for n > n that (cf. (15.30) of [11]) 

IK 1 

P[w"{U ni 8)>e]< — (F{t) - F(0))^F{25)) 2a - 1 (7.24) 

where the constant K' does not depend on 8. Given e > and n > 0, we can 
choose 8 so that 

(2K'/e i )(F(r) - F(0))(^(25)) 2 «" 1 < v , 

so by (7.24), we have (15.22) of [11] for n > no(8), and we can then follow [11] 
to get the desired tightness and convergence in D[0, r]. 

With the limiting Gaussian process denoted (Z(t),0 < t < r), we have from 
(7.23) and Theorem 3.2, along with the Skorohod coupling [38] and Fatou's 
lemma, that E[(Z(t) — Z(s)) 4 ] is bounded by Kg(t — s) 2a , so that by the 
Kolmogorov-Chentsov criterion (see e.g. [38]), the process G(-) admits a ver- 
sion with continuous sample paths. ■ 



Acknowledgements. Some of this work arose following discussions with 
Joseph Yukich. Travel for these discussions was partly funded by the London 
Mathematical Society. 

We thank the referee for pointing out a flaw in the presentation of Theorem 
2.1 in an earlier version of this paper. We thank Sana Louhichi and Xin Qi for 
providing preprint versions of their papers [18] and [54] repsectively. 



References 



[1] Andjcl, E. D. (1982) Invariant measures for the zero range processes. Ann. 
Probab. 10, 525-547. MR0659526 

[2] R. Atar, S. Athreya and M. Kang (2001) Ballistic deposition on a planar 
strip. Elec. Comm. Probab. 6, 31-38. MR1826165 

[3] Baddeley, A. and M0ller, J. (1989) Nearest-Neighbour Markov Point pro- 
cesses and random sets. Internat. Statist. Rev. 57, 89-121. 



M. D. Penrose/Particle systems: existence and limits 



31 



[4] Balazs, M., Rassoul-Agha, F., Seppalaincn. T. and Sethuraman, S. (2007) 
Existence of the zero range process and a deposition model with supcrlincar 
growth rates. Ann. Probab. 35 1201-1249. MR2330972 

[5] A.-L. Barabasi and H. E. Stanley (1995) Fractal Concepts in Surface 
Growth. Cambridge University Press. MR1600794 

[6] Baryshnikov, Y. and Yukich, J. E. Gaussian fields and random packing. J. 
Stat. Physics, 111, 443-463. MR1964280 

[7] Basis, V. Ya. Infinite-dimensional Markov processes with almost local inter- 
action of the components. Theory Probab. Appl. 21, 706-720. MR0426218 

[8] Basis, V. Ya. (1980) Stationarity and ergodicity of Markov interacting pro- 
cesses. Multicomponent random systems, pp. 37-58, Adv. Probab. Related 
Topics 6 Dckker, New York, 1980. MR0599532 

[9] J. D. Bcrnal (1959) A geometrical approach to the structure of liquids. 
Nature, 183 141-147. 
[10] J. D. Bernal and J. Mason (1960) Co-ordination of randomly packed 

spheres. Nature, 188, 910-911. 
[11] Billingslcy, P. (1968) Convergence of Probability Measures. Wiley, New 
York. MR0233396 

[12] Chen, M. F. (1992) From Markov Chains to Non-Equilibrium Particle Sys- 
tems. World Scientific, Singapore. MR2091955 
[13] H. S. Choi, J. Talbot, G. Tarjus and P. Viot (1993) First-layer formation in 

ballistic deposition of spherical particles: Kinetics and structure. J. Chem. 

Phys. 99, 9296-9303. 
[14] D. J. Cumberland and R.J. Crawford (1987) The Packing of Particles. 

Handbook of Powder Technology, 6, Elsevier, AMsterdam. 
[15] Davis, M. H. A. (1993) Markov Models and Optimization. Chapman and 

Hall, London. MR1283589 
[16] De Masi, A. and Presutti, E. (1992) Mathematical Methods for Hydrody- 

namic Limits. Springer, Berlin. 
[17] Dhar, D. (1999) The Abelian sandpilcs and related models. Phys. A 263, 

4-25. 

[18] Doukhan, P., Lang, G., Louhichi, S. and Ycart, B. (2005). A functional 
central limit theorem for interacting particle systems on transitive graphs. 
Markov Proc. Rel. Fields, to appear. Preprint available at ArXiv.math- 
ph/0509041. 

[19] R. Durrett (1988). Lecture Notes on Particle Systems and Percolation. 

Wadsworth & Brooks/Cole, Pacific Grove, CA. MR0940469 
[20] E. B. Dynkin (1965). Markov Processes, Volume L. Springer, Berlin. 
[21] J. W. Evans (1993) Random and cooperative adsorption. Rev. Modern 

Phys. 65, 1281-1329. 
[22] Ethier, S. N. and Kurtz, T. G. (1986) Markov Processes: Characterization 

and Convergence. Wiley, New York. MR0838085 
[23] Ferrari, P. A., Fernandez, R., and Garcia, N. L. (2002) Perfect simulation for 

interacting point processes, loss networks and Ising models. Stoch. Process. 

Applns. 102, 63-88. MR1934155 
[24] Fristedt, B. and Gray, L. (1997) A Modern Approach to Probability Theory. 



M. D. Penrose/Particle systems: existence and limits 



35 



Birkhauser, Boston. MR1422917 
[25] Garcia, N. L. and Kurtz, T. G. (2006) Spatial birth and death processes as 

solutions of stochastic equations. Alea 1, 288-303. MR2249658 
[26] Glotzl, E. (1981). Time-reversible and Gibbsian point processes I. Marko- 

vian and spatial birth and death processes on a general phase space. Math. 

Nachr. 102, 217-222. MR0642153 
[27] Glotzl, E. (1982). Time-reversible and Gibbsian point processes II. Marko- 

vian particle jump processes on a general phase space. Math. Nachr. 106, 

63-71. MR0675745 

[28] Griffcath, D. (1979). Additive and Cancellative Interacting Particle Sys- 
tems. Lecture notes in Math. 724, Springer, Berlin. MR0538077 

[29] A. Gervois, M. Lichtenberg. L. Ogcr, and E. Guyon (1989) Coordination 
number of disordered packings of identical spheres. J. Phys. A: Math, Gen. 
22 2119-2131. MR1004917 

[30] K. Gotoh and J. L. Finney (1974) Statistical geometrical approach to ran- 
dom packing density of equal spheres. Nature, 252, 202-205. 

[31] Harris, T. E. (1972). Nearest-neighbor Markov interaction processes on 
multidimensional lattices. Adv. Math. 9, 66-89. MR0307392 

[32] Holley, R. (1972) Markovian interaction processes with finite range inter- 
actions. Ann. Math. Stat. 43, 1961-1967. MR0373070 

[33] Holley, R. and Stroock, D. W. (1978) Nearest neighbour birth and death 
processes on the real line. Acta Math. 140, 103-154. MR0488380 

[34] Holley, R. ad Stroock, D. W. (1979) Central limit phenomena of various 
interacting systems. Ann. Math. 110, 333-393. MR0549491 

[35] Jacobsen, M. (2006) Point Process Theory and Applications: Marked Point 
and Piecewise Deterministic Processes. Birkhauser. MR2189574 

[36] W. S. Jodrey and E. M. Torey (1985) Computer simulation of close random 
packing of equal spheres. Phys. Review A, 32, 2347-2351. 

[37] R. Jullien and P. Meakin (1987) Simple three-dimensional models for bal- 
listic deposition and restructuring. Europhys. Lett. 4, 1385-1390. 

[38] Kallcnberg, O. (2002) Foundations of Modern Probability, 2nd ed., Springer, 
New York. MR1876169 

[39] Liggett, T. M. (1972) Existence theorems for infinite particle systems. 
Transactions Amer. Math. Soc. 165, 471-481. MR0309218 

[40] Liggett, T. M. (1985) Interacting Particle Systems. Springer, Berlin. 
MR0776231 

[41] Liggett, T. M. and Spitzer, F. (1981) Ergodic theorems for coupled random 
walks and other systems with locally interacting components. Z. Wahrsch. 
Verw. Gebiete 56 443-468. MR0621659 

[42] Maes, C. M., Rcdig, F. and Saada, E. (2002) The abelian sandpile model 
on an infinite tree. Ann. Probab. 30, 2081-2107. MR1944016 

[43] M0ller, J. and Waagepetersen, R. P. (2003). Statistical Inference and Sim- 
ulation for Spatial Point Processes. Chapman and Hall/CRC, Boca Raton. 
MR2004226 

[44] Onoda, G. and Liniger, E. (1990) Random loose packings of uniform spheres 
and the dilatancy onset. Phys. Rev. Letters, 64, 2727-2730. 



M. D. Penrose/Particle systems: existence and limits 



36 



[45] Penrose, M. D. (2001) Random parking, sequential adsorption, and the 
jamming limit. Comm. Math. Phys. 218, 153-176. MR1824203 

[46] Penrose, M. D. (2001) Limit theory for monolayer ballistic deposition in 
the continuum. J. Statist. Phys. 105, 561-583. MR1871657 

[47] Penrose, M. D. (2002) Limit theory for monotonic particle systems and 
sequential deposition. Stock Proc. Appln. 98 175-197. MR1887532 

[48] Penrose, M. D. (2005) Multivariate spatial central limit theorems with ap- 
plications to percolation and spatial graphs. Ann. Probab. 33, 1945-1991. 
MR2165584 

[49] Penrose, M.D. and Sudbury, A. (2005) Exact and approximate results for 

deposition and annihilation processes on graphs. Ann. Appl. Probab. 15, 

853-889 (2005). MR2114992 
[50] Penrose, M. D. and Yukich, J. E. (2002) Limit theory for random sequential 

packing and deposition. Ann. Appl. Probab. 12, 232-301. MR1890065 
[51] Penrose, M. D. and Yukich, J. E. (2001) Mathematics of random growing 

interfaces. J. Phys. A: Math. Gen. 34, 6239-6247. MR1862710 
[52] Preston, C. J. (1975) Spatial birth-and-death processes. Bull. Int. Statist. 

Inst. 46 (2) 371-391. MR0474532 
[53] Privman, V. (2000) Adhesion of Submicron Particles on Solid Surfaces, 

Colloids and Surfaces A, 165. Special Issue, edited by V. Privman. 
[54] Qi, X. (2007) Functional central limit theorem for spatial birth and death 

processes. Preprint, University of Wisconsin-Madison. 
[55] Scott, G. D. (1960) Packing of spheres. Nature, 188, 908-909. 
[56] B. Senger, J.-C. Voegel, P. Schaaf (2000) Irreversible adsorption of colloidal 

particles on solid substrates. Colloids and Surfaces A 165, 255-285. 
[57] Seppalalncn, T. (2000) Strong law of large numbers for the interface in 

ballistic deposition. Ann. Inst. H. Poincare Probab. Statist. 36, 691-736. 

MR1797390 

[58] Sidoravicius, V. and Vares, M. E. (1995) Ergodicity of Spitzer's renewal 
model. Stochastic Process. Appl. 55 119-130. MR1312152 

[59] H. Solomon (1967) Random Packing Density. Proc. Fifth Berkeley Symp. 
on Prob. and Stat. 3, 119-134, Univ. California Press. MR0256531 

[60] Spohn, H. (1991) Large scale dynamics of interacting particles. Springer, 
Berlin 

[61] T. Vicsek (1989) Fractal Growth Phenomena. World Scientific, Singapore. 
MR1020762 

[62] Williams, D. (1991) Probability with Martingales. Cambridge University 

Press, Cambridge. MR1 155402 
[63] Zong, C. (1999) Sphere Packings. Springer, New York. MR1707318 



